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We observe interlaced square vortex lattices in rotating dilute-gas spinor Bose-Einstein conden- 
sates (BEC). After preparing a hexagonal vortex lattice in a one- component BEC in an internal 
atomic state |1), we coherently transfer a fraction of the superfluid to a different state |2). The 
subsequent evolution of this pseudo-spin- 1/2 superfluid towards a state of offset square lattices in- 
volves an intriguing interplay of phase-separation and -mixing dynamics, both macroscopically and 
on the length scale of the vortex cores, and a stage of vortex turbulence. The stability of the square 
structure is proved by its response to applied shear perturbations. An interference technique shows 
the spatial offset between the two vortex lattices. Vortex cores in either component are filled by 
fluid of the other component, such that the spin- 1/2 order parameter forms a Skyrmion lattice. 

PACS numbers: 03.75.Lm,03.75.Mn,67.90.+z,73.43.-f,71.70.Di,67.40.Vs,32.80.Pj 



Some of the most intriguing phenomena of superfluid- 
ity are revealed when a superfluid is set rotating. De- 
pending on the complexity of the order parameter a su- 
perfluid can carry angular momentum in different ways. 
A single-component superfluid exhibits quantized vor- 
tices that, in a dilute-gas single-component BEC, orga- 
nize into a regular hexagonal Abrikosov lattice jl| [see Fig. 
1(a)], whose static properties and long- wavelength exci- 
tations have been examined [2, lJlIIlJ- Such lattices are 
also formed in Type II superconductors in a magnetic 
field 1] and in superfluid 4 He Q. Multi-component su- 
perfluids, with their pseudo-spin order parameters, carry 
angular momentum in spin textures (Skyrmions) 0, , 
which may again form lattices 0,0, El ? such as in the 
superfluid phases of liquid 3 He |l2j and in certain regimes 
of two-dimensional electron systems . 

In this work we study rotating two-component super- 
fluids in a magnetically trapped dilute-gas BEC Q de- 
scribed by a pseudo-spin- 1/2 order parameter. State |1) 
represents spin-down and state |2) represents spin-up, 
and the azimuthal component of the spin is fixed by the 
relative phase $ between the state |1) and |2) wavefunc- 
tions. Both the large-scale density profiles and the vor- 
tex lattice structure depend on the nature of inter- and 
intracomponent interactions. The interaction energy, de- 
termined by three scalar coupling constants gn, #22 and 
gi2 and the population densities n\ and n2, is given by 
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For attractive intercomponent interaction (gi 2 < 0) it 
is energetically favorable if the spatial profiles of both 
components coincide, and a hexagonal vortex lattice is 
expected, as in the one-component case. For repulsive in- 
tercomponent interaction (gi2 > 0), realized in our 87 Rb 
system, it is energetically favorable for the two compo- 
nents to reduce spatial overlap. In a static BEC this may 
be achieved by macroscopic phase separation, whereas in 
rotating two-component BEC the presence of the vortex 




FIG. 1: (a) Hexagonal vortex lattice in a one-component 
BEC. (b) Square lattice, viewed in the |1) component of a two- 
component BEC. (c) Reciprocal (k— ) space of the hexagonal 
lattice, obtained by 2D Fourier transform of (a), showing 6 
peaks spaced by 60°. (d) k— space of the square lattice. The 
reciprocal lattice vectors ki and k^ enclose an angle (£12 = 90°. 



lattice allows for more subtle separation effects, such as 
interlacing the lattices. Vortices in either component may 
be filled by the other component, and the pseudo-spin- 
1/2 order parameter forms a Skyrmion lattice. Studies 
by Mueller and Ho 9] and Kasamatsu et al. [l(| predict, 
among other structures, interlaced square lattices over a 
wide range of parameters. Presumably this is because in 
a square lattice, unlike in a triangular lattice, each vor- 
tex in component |1) can have all its nearest-neighbor 
vortices be in component |2), and vice versa |l5|- 

Under suitable conditions we indeed observe a spatial 
separation of vortices in the two components, followed by 
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the formation of an ordered square lattice structure, as 
shown in Fig. 1(b). Its formation and decay dynamics, 
as well as its static and dynamic properties are examined 
in this work. 

Our 87 Rb two-state system Q consists of two mag- 
netically trappable hyperfine-Zeeman levels of the 87 Rb 
atom - \F = l,mF = —1), henceforth called |1), and 
\F = 2,mf = 1), henceforth called |2). Spatial over- 
lap of the two states is realized in a harmonic, axi- 
ally symmetric magnetic trap with oscillation frequencies 
{cu p ,u; z } = 27r{7.7, 4.9} Hz. An electromagnetic coupling 
between the two states can be achieved via a two-photon 
transition, involving a microwave photon at ~ 6.833 GHz 
and a radio frequency photon at ~ 1 MHz. A short-pulse 
application of this coupling drive yields a desired, spa- 
tially uniform amount of population transfer between the 
two states, instantaneous with respect to the external dy- 
namics of the two states. In the pseudo-spin- 1/2 picture 
the effect of the coupling drive is to cause spin-rotations 
[l^ . Inelastic atomic collisions limit the lifetime of the |2) 
population to a few seconds. The decay does not cause 
|2) atoms to convert back to the |1) state but rather to 
leave the trap altogether. 

To study the rotational properties of this two- 
superfluid system, we initially create regular hexagonal 
vortex lattices in BECs in state |1), as described in ear- 
lier work 0, 0, • Here we start with near-pure conden- 
sates containing (3.5 — 4) x 10 6 atoms, rotating at a rate 
ft w 0.75 x uo p about the z-axis. A "transfer pulse" of 
the coupling drive then transfers a fraction of the popu- 
lation into state |2) , in this work 80 — 85 % as discussed 
below. After a variable wait time we take two images 
of the system. A nondestructive phase-contrast image is 
taken either of one component alone or of both compo- 
nents simultaneously |l4| , along an axis perpendicular to 
the rotation axis ("side view") while the system is still 
trapped. In order to optically resolve the vortex struc- 
tures we expand either the |1) or |2) component of the 
condensate by a factor of 9, to a diameter of ~ 600 /im, 
before a second, destructive image is taken along the ro- 
tation axis ("top view"). The other component is re- 
moved at the beginning of the expansion. The details of 
expansion and imaging have been described in Ref. [3. 

In Fig. 2 we analyze the formation and decay dynam- 
ics of the square lattice structure. Figs. 2(a) and (b) 
show time sequences of top- view images of the expanded 
|1) and |2) states, taken in different experimental runs. 
The time evolution after the transfer pulse involves sev- 
eral stages. For the first ~ 0.1 — 0.25 s surprisingly little 
dynamics is visible and certainly no structural transi- 
tion in the vortex lattice is seen in either component. 
From ~ 0.25 — 2 s a turbulent stage evolves in both com- 
ponents in which vortex visibility degrades significantly, 
shown in Fig. 2[a(ii)] and Fig. 2[b(ii)]. As we will show, 
this turbulence is directly linked with the transition from 
overlapping hexagonal vortex lattices to interlaced square 
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FIG. 2: (a) Time sequence of images of state |1), after ~ 80% 
population transfer to |2), showing evolution from a hexagonal 
lattice over a turbulent stage to a square structure and back 
to a hexagonal lattice, (b) Also state |2) forms a square lattice 
(iv) before its decay, (c) Detailed time evolution in reciprocal 
space. Intensity in an annulus along the ip coordinate [defined 
in Fig. 1(c)] is shown on the ordinate. The initial 6-peak 
structure of the hexagonal lattice vanishes quickly because of 
turbulence. From 3 — 5.5 s square lattices give rise to a 4- 
peak structure. Around 6 — 9 s a transition back to hexagonal 
lattices occurs, (d) Two-color side view images of the initial 
turbulent evolution. State |1) (|2)) appears bright (dark) on 
gray background. The fine filament structures in (ii-v) are 
due to mutual filling of vortex cores. 



lattices. From 2 — 3 s square domains emerge from the 
turbulent state, and defects propagate out of the lattice 
{Fig. 2[a(iii)]}. From 3 — 5.5s stable square lattices 
are observed in both components {Fig. 2[a(iv)] and Fig. 
2[b(iv)]}. At this stage, around 4s, despite the large 
(80 — 85%) initial population transfer to state |2), the 
number of |2) atoms has decreased to only 1.5 x 10 5 , 
while state |1) contains (5 — 7) x 10 5 atoms. As the |2) 
state population continues to decay, the vortex lattice 
planes bend {Fig. 2[a(v)]} and a transition back to a 
hexagonal lattice in state |1) takes place {Fig. 2[a(vi)]}. 
During the transition from square lattices to a hexagonal 
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lattice no turbulence occurs. 

A more quantitative analysis of these dynamics is pos- 
sible in reciprocal space. Figure 2(c) shows the time evo- 
lution of the intensity within an annulus in reciprocal 
space, defined in Fig. 1 (a), which contains the recip- 
rocal lattice peaks. Initially six peaks are visible, sepa- 
rated by 60°, forming the reciprocal lattice of a hexagonal 
vortex lattice. Because of turbulence these peaks vanish 
between ~ 200 ms and 2 s. After 2 — 3 s a four-peak struc- 
ture appears, that is stable for a period of ~ 2.5 s. The 
observed ratio of reciprocal lattice vector lengths [defined 
in Fig. 1(d)] k 1 /k 2 = 0.98(2) and the angle cp 12 = 95(3)° 
between k\ and k 2l clearly identify a square lattice. At 
~ 5.5 s the appearance of two additional peaks signals the 
onset of a transition back to the hexagonal state, which 
is completed by ~ 9 s. 

To examine the origin of the initial turbulence, we 
show in Fig. 2(d) a time sequence of in-trap side-view 
images, where the |1) (|2)) state appears bright (dark). 
As visible in Fig. 2[d(i)], a macroscopic component sep- 
aration from the initially homogeneous two-component 
superposition to a ball-shell structure takes place within 
50 — 100 ms 14]. During this period of dramatic axial 
separation, both the individual vortices and the overall 
vortex lattice remain remarkably quiet, as viewed along 
the rotation axis {see Fig. 2[a(i)] and Fig. 2[b(i)]}. 
Around 600 ms, fine filament structures appear at the in- 
tercomponent boundary {Fig. 2[d(ii)]}, coincident with 
the full development of vortex turbulence seen in top- 
view images {Fig. 2[a(ii)]}. The filament structures dis- 
tort and fill out the whole BEC as the vortex turbulence 
peaks at around Is {Fig. 2[d(iii)]}, and straighten up 
at around 2 s {Fig. 2[d(iv)]} coincident with the restora- 
tion of vortex visibility in top view images. Subsequently 
the filaments become less visible as state |2) decays {Fig. 
2[d(v)]}. We interpret these structures as vortex cores 
in either component being filled by fluid of the other 
component, forming a Skyrmion lattice Filled vortex 
cores grow in size above the resolution limit of our side 
view images and become observable, while empty vortices 
in single-component BECs in equilibrium are well below 
this resolution limit and are not observed in-trap [l7^ . 
The initial macroscopic, predominantly axial phase sep- 
aration has thus evolved into a microscopic separation 
of two interlaced vortex lattices. It is this microscopic 
separation which gives rise to the observed turbulence. 

To study the stability of the square lattice , we excite 
shear perturbations (Tkachenko modes) |l8|, LL9| in the 
square lattice by focusing a resonant laser beam onto the 
center of the condensate 3]. As shown in Fig. 3, the 
perturbation relaxes to the equilibrium square configu- 
ration within 500 — 800 ms after excitation. However, in 
contrast to single-component triangular lattices |3|, no 
clear oscillation is observed. Two effects may contribute: 
The excitation may be overdamped because of ill-defined 
boundary conditions caused by imperfections in the outer 
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FIG. 3: (a) Tkachenko excitation in the square lattice, ob- 
served 250 ms after beginning of excitation, and (b) after re- 
laxation back to a square lattice structure (at 850 ms). 




FIG. 4: (a) Vortex lattice interference, (i) and (iii): Simula- 
tions for different values of the relative phase <I> between the 
|1) and |2) state wavef unctions, (ii) and (iv): Experimental 
results, (b) Time evolution of the interference patterns shows 
no vortex offset before the turbulence (i) but patterns char- 
acteristic of offset vortices afterwards [b(iii)],[a(ii)],[a(iv)]. 

region of the square lattice, where the |2) state popula- 
tion has already decayed. The oscillation may also be 
masked by random lattice excitations of comparable am- 
plitude that cannot be completely removed in the short 
time between formation and decay of the square lattice. 
However, the return of the lattice to its square config- 
uration is sufficient to demonstrate the stability of the 
square structure in the two-component system. 

So far we have presented evidence that both compo- 
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nents separately form regular and stable square lattices. 
Filament structures indirectly indicate a microscopic- 
scale spatial separation of vortices. In the following, we 
employ an interference technique between the two su- 
perfluids to more concisely address two questions: Are 
the vortex lattices really offset from each other? Do the 
vortices really only separate from each other during the 
turbulent stage, after a surprising delay of 200 ms? Fig- 
ures 4[a(i)] and (hi) show results of a simple simulation of 
interference between two square vortex lattice wavefunc- 
tions, spatially offset by 1/2 (a + b), where a, b are the 
two basis vectors of the square lattice. Two values (<£ = 
and = 45°) of the relative quantum phase $ between 
the two wavefunctions 20] are considered. A population 
ratio of A^| 2 > : A^|!> = 20 : 80 was assumed, as in the ex- 
periment at 4 s. Evidently two very dissimilar patterns 
are seen, consisting either (<£ ~ 0) of dark patches in be- 
tween two vortices, surrounded by a bright square struc- 
ture, or (<I> ~ 45°) of staircase-like vortex arrangements 
separated by bright stripes. We checked that the quali- 
tative appearance of such simulations changes drastically 
when altering the spatial offset from 1/2 (a + b). 

To observe such structures experimentally, we apply 
a 7r/2 "interference" pulse to the two-component sys- 
tem just before expansion for the top view image. This 
pulse results in a phase-coherent transfer of population 
between the two states, thus creating interference. If the 
7r/2 pulse is applied under conditions when regular square 
lattices are expected in both components, the observed 
images, Fig. 4[a(ii)] and (iv), agree qualitatively with the 
simulations, clearly demonstrating a spatial offset close 
to l/2(a + b) between the vortex lattices. 

The qualitative agreement between experimental re- 
sults and simulation indicates that vortex lattice interfer- 
ence may be employed to examine the separation process 
of the vortex lattices. When varying the wait time be- 
tween the transfer pulse and the tt/2 interference pulse, 
we observe the following qualitative features [Fig. 4(b)]: 
Until turbulence occurs (after 200 — 300 ms) there are 
no indications for separation of the vortex lattices. Fig- 
ure 4.(b)(i) shows a hexagonal lattice with good vortex 
contrast, very similar to images of each single compo- 
nent at this time. This confirms that during the initial 
200 — 300 ms of dramatic macroscopic component sepa- 
ration, vortices in the two components continue to form 
identical and overlapped lattice structures, when viewed 
along the rotation axis. Only after this surprisingly long 
delay, and after the turbulent stage do interference pat- 
terns show a grouping characteristic of offset vortices 
{Fig. 4[b(iii)]}. This observation is further proof of the 
direct link between the turbulent period and the micro- 
scopic separation of the vortex lattices. 

In conclusion, we observe a new vortex lattice struc- 
ture in rotating two-component BEC. Each component 
carries a square vortex lattice, and the lattices are inter- 
laced. Vortices in both components are filled by fluid of 



the other component, forming a Skyrmion lattice. This 
structure is stable, as evidenced by relaxation of applied 
shear excitations back to the square structure. An in- 
triguing turbulent stage accompanies the transition from 
overlapped hexagonal lattices to interlaced square lat- 
tices. Vortex lattice interference is used to study the 
offset dynamics of the two lattices. 
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